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Abstract: Analogue spacetimes, (and more boldly, analogue models both of and for 
gravity), have attracted significant and increasing attention over the last decade and a 
half. Perhaps the most straightforward physical example, which serves as a template 
for most of the others, is Bill Unruh's model for a dumb hole, (mute black hole, acous- 
tic black hole), wherein sound is dragged along by a moving fluid — and can even 
be trapped behind an acoustic horizon. This and related analogue models for curved 
spacetimes are useful in many ways: Analogue spacetimes provide general relativists 
with extremely concrete physical models to help focus their thinking, and conversely 
the techniques of curved spacetime can sometimes help improve our understanding of 
condensed matter and/or optical systems by providing an unexpected and counter- 
vailing viewpoint. In this chapter, I shall provide a few simple examples of analogue 
spacetimes as general background for the rest of the contributions. 
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1 Introduction 

While the pre-history of analogue spacetimes is quite long and convoluted, with optics- 
based contributions dating as far back as the Gordon metric of 1923 [1], significant 
attention from within the general relativity community dates back to Bill Unruh's 
PRL concerning acoustic black holes (dumb holes) published in 1981 [2]. Even then, it 
is fair to say that the investigation of analogue spacetimes did not become mainstream 
until the late 1990's. (See the recently updated Living Review article on "Analogue 
gravity" for a summary of the historical context [3].) 

In all of the analogue spacetimes, the key idea is to take some sort of "excitation", 
travelling on some sort of "background", and analyze its propagation in terms of the 
tools and methods of differential geometry. The first crucial technical distinction one 
has to make is between "rays" and "waves" . 
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• The rays of ray optics (geometrical optics), ray acoustics (geometrical acoustics), 
or indeed any more general ray-like phenomenon, are only concerned with the 
"light cones" , "sound cones" , or more generally the purely geometrical "propaga- 
tion cones" defined by the ray propagation speed relative to the appropriate back- 
ground. Physically, in this approximation one should think photons/phonons/ 
quasi-particles following some well-localized trajectory, rather than the more dif- 
fuse notion of a wave. Mathematically, we will soon see that it is appropriate to 
construct some metric g a b, and some tangent vector k a to the particle trajectory, 
such that: 



Here indices such as a, b, c, ... , take on values in {0, 1, 2, 3}, corresponding to 
both time and space, whereas indices such as i, j, k, . . . , will take on values in 
{1, 2, 3}, corresponding to space only. Of course we could multiply the metric by 
any scalar quantity without affecting this equation; this is known as a conformal 
transformation of the metric. (So distances change but angles are unaffected.) 
In the language of differential geometry, ray phenomena are sensitive only to a 
conformal class of Lorentzian geometries. 

• In contrast, for waves one needs to write down some PDE — some sort of wave 
equation. For example, for a scalar excitation \I/ one needs to construct a wave 
equation in terms of a d'Alembertian [2-7]: 



This d'Alembertian, (and in fact very many of the different possible types of wave 
equation), depends on all the components of the metric g a b, not just the conformal 
class. (And conformal wave equations, of which the Maxwell electromagnetic wave 
equations are the most common, have their own somewhat different issues.) 

In short, depending on exactly what one is trying to accomplish, one may sometimes 
be able to get away with ignoring an overall multiplicative conformal factor — but for 
other applications knowledge of the conformal factor is utterly essential. 

What I shall now do is to present some elementary examples — and a few not so 
elementary implications — that will hopefully serve as a pedagogical introduction to 
the more specific physics problems addressed in the other contributions to this volume. 



gab k a k b = 0. 
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2 Optics: the Gordon metric and its generalizations 



The original Gordon metric [1] of 1923 was limited to ray optics in a medium with 
a position-independent refractive index, and with some position-independent velocity. 
Let 



Vab 



-1 








Sij 



(2-1) 



denote (as usual) the special relativistic Minkowski metric, and correspondingly set 
the zeroth coordinate to x° = t = ct p h y sicai- Denote the refractive index by n and the 
4-velocity of the medium by V a = 7(1; /3n). Then we have V a = j(— 1; /3n). Now 
define 

9ab = (Vab + V a V b ) ^ = Tjab + 



n 



1 - 



V a V b . 



(2.2) 



In the rest frame of the medium V a — > (1; 0) and 



9ab -»■ 



-1/n 2 
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(2.3) 



Therefore in this rest frame the null cones of the medium are exactly what we want: 

dt 2 



= ds" 



+ dx 



rr 



dx 
dt 



1 

n 



(2.4) 



But more generally, for non-zero velocity, 7^ 0, the metric g ab provides a perfectly 
good special relativistic model for the light cones in a homogeneous moving medium. 
Let us agree to raise and lower the indices on the 4-velocity V using the Minkowski 
metric rj, then the contravariant Gordon metric is 



9 



ab 



(r] ab + V a V b ) - n 2 V a V b 



(2.5) 



The first and most obvious generalization is to note that one can easily make 
the refractive index and 4-velocity both space and time dependent. (Physically, this 
will certainly work as long as the wavelength and period of the light wavicle is short 
compared to the spatial and temporal scales over which changes of the background 
refractive index and 4-velocity are taking place.) 

A second generalization is to note that from the point of view of ray optics one 
might as well take 



9ab 



(Vab + VaVb) 



VaVb 



n- 



g ab = Q- 2 [( v ab + V a V b ) - n 2 V a V b ] . (2.6) 
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The conformal factor will simply drop out when determining the light cones. With 
the quantities Q(x), n(x), and V(x) all being space and/or time dependent, this is the 
most general (but still physically natural) form of the (special relativity based) Gordon 
metric one can write. 

One can certainly calculate the Einstein tensor for this optical metric, but there 
is a priori no really compelling reason to do so — there is a priori no good reason 
to attempt to enforce the Einstein equations for this optical metric, the physics is 
just completely different. (That being said, if one merely views this as an ansatz 
for interesting metrics to look at, then many of the standard spacetimes of general 
relativity can certainly be put into this form. For example, the Schwarzschild and 
Reissner-Nordstrom spacetimes, and the FLRW cosmologies, can certainly be put in 
this form [8].) 

Example: Let us take Q = n and write 

g ab = n 2 { Vab + V a V b ) - V a V b . (2.7) 

This corresponds to 

ds 2 = -dt 2 + n 2 ||dx 2 ||. (2.8) 
Now pick the specific refractive- index profile 

n =T^T-2> ( 2 - 9 ) 
1 + r z /a z 

so 



,, \dr 2 + r 2 {d6 2 + sm 2 6 d(j) 2 }] 
ds 2 = -dt 2 + -5i — ^ - 9N9 (2.10) 



(l + r 2 /a 2 ) 2 

A theoretical cosmologist should recognize this as the Einstein static universe in isotropic 
coordinates [9-11]. A theoretician working in optics should recognize this as the 
Maxwell fish-eye lens [12-14]. This is simply the first of many cross-connections be- 
tween optics and general relativity. This becomes (or should become) a two-way street 
for information exchange. 

The Maxwell fish-eye above is an example of a Liineburg lens [15], and has now 
become the canonical example which helped initiate much of the recently developed 
field of "transformation optics" [16-19]. In particular, if one looks at this from the 
perspective of a theoretical cosmologist then the prefect focussing properties are utterly 
trivial — after all, the spatial slices of the Einstein static universe are just the hyper- 
sphere S 3 in suitable coordinates — the geodesies are obviously just great circles, which 
by symmetry must meet at the antipodes of the emission event, and so perfect focussing 
in the ray optics approximation is trivial. 
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Limitations: Perhaps the greatest limitation of the Gordon metric is its inability (in 
its original 1923 formulation) to deal with wave properties of light. There is a rather 
non-trivial generalization to the full Maxwell equations [3], but for technical reasons 
the generalization requires the very specific constraint 



[magnetic permittivity] oc [electric permeability]. 



(2.11) 



For ordinary physical media this constraint is somewhat unphysical [3], but there is hope 
that suitably designed metamaterials [20] may be designed to at least approximately 
satisfy this constraint. 

Foreground-background version: So far, the Gordon metric has been based on a 
optical medium in Minkowski space described by the flat metric r] ab . But now suppose 
we have a non-trivial background metric f a b arising from standard general relativity, 
and place a flowing optical medium on top of that. It now becomes interesting to 
consider the generalized Gordon metric 



9ab 



Q 2 



(fab + V a V b ) 



VnVh 



n- 



g ab = Q- 2 [(f ab + V a V b ) - n 2 V a V b ] . (2.12) 



We now have the possibility of a non-trivial general relativistic background f a b(x), a 
position-dependent refractive index n(x), a position-dependent 4- velocity V a (x), and 
a position-dependent conformal factor Q(x). Note that the 4- velocity V a (x) now has 
to be a timelike unit vector with respect to the background metric f a b(x), and the 
indices on V are raised and lowered using /. In particular, note g a bV a V b — —Q 2 /n 2 
and g ab V a Vb = —n 2 /Q 2 , so it makes sense to define 

Q 



ya 



n 



ya 



and 



Vn 
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v a . 



(2.13) 



Then V is a timelike unit vector with respect to g, and its indices should be raised and 
lowered using g. Then we can adapt the generalized Gordon metric to also write the 
background / in terms of the foreground g as: 



ah 



n 



-2 



(gab + V a V b ) - n 2 V a V b 



f 



ab 



Q 2 



(g ab + V a V b ) - 



yayb 



ir 



(2.14) 



For a relativist the generalized Gordon metric provides an interesting ansatz for a 
potentially intriguing class of spacetimes to consider. From a theoretical optics per- 
spective, one might view this procedure as an extremely general way of "composing" 
and/or "inverting" the transformations of transformation optics — for example, one 
might first design some metamaterial [20] to generate the background f a b(%), and then 
impose some flowing optical medium on top of that. Various interesting possibilities 
come to mind. 
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3 Non-relativist ic acoustics: the Unruh metric 



Bill Unruh's 1981 PRL article [2], and much of the follow up work [3-5], was explicitly 
and intrinsically based on non-relativistic acoustics. Let us explore the basic features 
of this particular model. 

Geometric acoustics: From the acoustic ray perspective the derivation is trivial: 
Let c s be the speed of sound, and let v be the velocity of the fluid. Then sound rays 
(phonon trajectories) satisfy [3-5] 

||dx- v dt|| = c s dt. (3.1) 

Let us, already anticipating the possibility of an arbitrary conformal factor, define 

ds 2 = tt 2 {-c s dt 2 + (dx - v dt) 2 } = tt 2 {-(c 2 - v 2 ) - 2 v ■ dx dt + ||dx|| 2 } . (3.2) 

(The zeroth coordinate is now most naturally chosen to simply be x° = t = £ p h yS icai, 
without any explicit factor of c. The speed of sound c s has the dimensions of a physical 
velocity.) Then the sound-ray condition is completely equivalent geometrically to the 
null-cone condition ds 2 = 0. In terms of a 4 x 4 matrix this is equivalent to defining 
the metric tensor [3-5] 



9ab = ^ 2 

The corresponding inverse metric is 







-(c 2 -v 2 ) 









.ah 
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-l/c 2 


-vi/c 2 


-v^c 2 





(3.3) 



(3.4) 



It should be emphasized that in this situation the velocity v and speed of sound c s will 
be inter-related in some (often quite complicated) manner — the background fluid flow 
must satisfy the Euler equation and the continuity equation [2-5]. 

It should again further be emphasized (forcefully) that while one can certainly 
calculate the Einstein tensor for this acoustic metric, there is a priori no really com- 
pelling reason to do so — there is a priori no good reason to attempt to enforce the 
Einstein equations for this acoustic metric, the physics is just completely different. 
That being said, if one again views this as an ansatz for interesting metrics to look 
at, then many of the standard spacetimes of general relativity (but certainly not all 
interesting spacetimes) can be put into this form. (For instance the Schwarzschild 
and Reissner-Nordstrom spacetimes can be put into this form by going to Painleve- 
Gullstrand coordinates, but the Kerr and Kerr-Newman spacetimes cannot be put in 
this form [3, 21].) 



- 6 - 



To further develop the discussion, let us now introduce quantities 



Q 



ab 















V a = (1; v') = (1; v). 



(3.5) 



Here the 4-velocity V a is normalized non-relativistically — with the time component 
being unity. Then for the inverse metric 
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ab 



Q ab - 



yayb 



(3.6) 



But what about the covariant metric q nh 7 Let us now define 



Qab 












Sij 



(3.7) 



Then Q b ab is the Moore-Penrose pseudo- inverse of Q ab , and the object 



P\ = Q ac Q 



do 















(3.8) 



is a projection operator onto spatial slices. Let us now furthermore define the quantities 
V* = Q b ab V b = (0; v*) = (0; v), while T a = (1; 0). The best we can do for the covariant 
metric g ao is to now write the somewhat clumsy expression: 



9ab = [QL 



T a T b -T a Vt-V*T b }. 



(3.9) 

2 



In view of the fact that, with these definitions, one has T a V a = 1 and V a = v 
while Q ab T b = and Q ab V b = V a — T a , it is easy to verify that (as required) 



9 9bc 



S a 



(3.10) 



As we shall soon see, relativistic acoustics is in some sense actually somewhat simpler 
than the non-relativistic case. 

Wave acoustics: If one goes beyond ray acoustics, then the parameter Q is no longer 
arbitrary. One does have to make some additional (and rather stringent) technical 
assumptions — barotropic, irrotational, and inviscd (zero viscosity) flow [3-5]. Under 
those assumptions, by linearizing the Euler equation and continuity equation, after a 
little work one ultimately obtains a wave equation (a curved-spacetime d'Alembertian 
equation) for perturbations of the velocity potential specified in terms of the density 
of the fluid and the speed of sound — specifically one has = ^/p/c s in 3 space 
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dimensions, Q = p/c s in 2 space dimensions, and technical problems arise in 1 space 
dimension. Generally, in d space dimensions, Vt = (p/c s ) 1 ^ d ~ 1 \ (See for example 
reference [3].) 

A specific feature of physical (wave) acoustics, not probed in the geometrical acous- 
tics limit, is the behaviour of quasi-normal modes [22, 23]. Furthermore, if the flow is 
not irrotational, so one is dealing with both background vorticity and vorticity-bearing 
perturbations, then a considerably more complicated system of wave equations can be 
written down [24], but this system of PDEs has nowhere near as clean a geometrical 
interpretation as the irrotational case. 

4 Horizons and ergo-surfaces in non-relativistic acoustics 

One of the very nice features of non-relativistic acoustics is that it is very simple and 
straightforward to define horizons and ergo-surfaces [3-5]. To define these concepts, 
it is sufficient to work in the geometric acoustics limit; wave acoustics adds additional 
constraints not directly needed to define horizons and ergo-surfaces. Consider for sim- 
plicity a stationary (time independent) configuration. 

• Ergo-surfaces are defined by the condition ||v|| = c s . 

• Horizons are surfaces, located for definiteness at /(x) = 0, that are defined by 
the 3-dimensional spatial condition V/ • v = c s ||V/||. 

So the ergo-surface bounds the region where one cannot stand still without generating 
a sonic boom, and corresponds to Mach one, (M = v/c s = 1). In contrast, on a horizon 
the normal component of the fluid velocity equals the speed of sound, thereby either 
trapping or anti-trapping the acoustic excitations. 

Stationary versus static: In general relativity the words "stationary" and "static" 
have precise technical meanings that may not be obvious to non-experts. So a few 
words of explanation are called for: 

• Stationary: For all practical purposes this means "time independent". More 
precisely, mathematically there is a Killing vector (a symmetry of the system) 
which is timelike at spatial infinity. Physically there is a class of natural time 
coordinates (not quite unique) in which the metric is time- in dependent. In this 
coordinate system the Killing vector is naturally associated with invariance under 
time translations t — > t + C. 
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• Static: For all practical purposes this means "time independent and non-rotating" . 
More precisely, mathematically there is a Killing vector (a symmetry of the sys- 
tem) which is both timelike at spatial infinity and "hypersurface orthogonal", 
meaning there exist functions £(x) and r(x) such that K a = £(rr) g ab dbr(x). 
Physically there is then a unique natural time coordinate, (in fact r, which is un- 
fortunately not necessarily "laboratory time"), in which the metric is both time- 
independent and block-diagonal. That is, with vanishing time-space components 
gu = 0, in these coordinates the metric block diagonalizes into (time) © (space). 
The existence of a coordinate system with vanishing time-space metric compo- 
nents is sufficient in general relativity to imply zero angular momentum for the 
spacetime, and absence of "frame dragging" , hence the sobriquet "non-rotating" . 

A word of warning: Just because one can always choose a coordinate system to block 
diagonalize a static spacetime does not mean this is always a good idea. Coordinates 
in which static spacetimes are block diagonal will break down at any horizon that 
might be present in the spacetime. (For instance, Schwarzschild geometry in the usual 
coordinates.) Permitting coordinates for static spacetimes which retain the manifest 
time independence, but do not explicitly force block diagonalization of the metric, 
has significant technical and physical advantages. For one thing, this is the most 
natural situation when one works with "laboratory time" and a time independent fluid 
flow. For another thing, once one allows off-diagonal elements for the metric one can 
easily construct "horizon penetrating" coordinates, which are well defined both at and 
across the horizon. (For instance, Schwarzschild geometry in Painleve-Gullstrand or 
Eddington-Finklestein coordinates.) In particular, the acoustic metric as given above 
(in terms of laboratory time, speed of sound, and fluid velocities) is automatically in 
horizon-penetrating form, all the components of both the metric g a b and its inverse g ab 
remain finite as one crosses the horizon. Let us now see how these ideas are used in 
practice. 

Static configurations: Suppose the background flow satisfies the integrability con- 
straint 

^^=V$, (4.1) 
cj — v 2 

and then consider the new time coordinate r = t + $. (Here t is explicitly laboratory 
time, while r is constructed for mathematical convenience rather than for direct physical 
purposes.) Note that this integrability condition implies (but is stronger than) the 
vanishing of local helicity 

h = v • (V x v) = 0. (4.2) 
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In terms of this new time coordinate 



ds 2 



^{-{cl-v 2 ) dr 2 



ViVj 



(4.3) 



The geometry is now in this form block-diagonal so it is manifestly static, not just 
stationary. (And so the ergo-surfaces and horizons will automatically coincide.) The 
time translation Killing vector is 

K a = (1; 0), so K a = -n 2 (c 2 s - v 2 ; 0). (4.4) 

To explicitly verify that this is hypersurface orthogonal in the sense defined above, note 



K a = -Q 2 {c 2 s - v 2 ) d a r, so K a = ■ 
The norm of this Killing vector is given by 

K a K a = -n 2 (c 2 -v 2 ) = -n 2 c 2 



(4.5) 



(4.6) 



Using very standard techniques, the surface gravity is then calculable in terms of the 
gradient of this norm [3]. It is a standard result that for a Killing horizon the over- 
all conformal factor drops out of the calculation [25]. Generalizing Unruh's original 
calculation [2], which corresponds to c s being constant, one finds [3-5] 



9h 



n • Vict 



v 2 ) 



H 



n ■ V(c„ 



H 



d(c s - v) 



dn 



(4.7) 



H 



which can also be compactly written in terms of the Mach number M = v/c s as 













DM 


9h = c 2 H 


n ■ V(v/c s ) 




n- VM 


, = 4 


dn 



(4- 



H 



We emphasise that this already works in the geometric acoustics framework, and that 
there is no need to make the more restrictive assumptions corresponding to wave acous- 
tics that were made in references [3] and [5]. If the integrability condition is not satisfied 
one must be a little more devious. 

Stationary but non-static configurations: If the acoustic horizon is stationary 
but not static there may or may not be additional symmetries, (in addition to the 
assumed time independence), so in particular the horizon may or may not be a Killing 
horizon. (A horizon is said to be a Killing horizon if and only if there exists some 
Killing vector such that the location of the horizon coincides with the vanishing of the 
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norm of that Killing vector. So Killing horizons automatically satisfy nice symmetry 
properties.) For a Killing horizon the calculation of surface gravity is still relatively 
straightforward, for non-Killing horizons the situation is far more complex. 

Note that in full generality, on the horizon we have (V/ • v) 2 = c 2 ||V/|| 2 , which 
we can rewrite in 3-dimensional form as g 1 ^ dif djf = 0. Since the configuration, and 
location of the horizon, is time independent this statement can be bootstrapped to 3+1 
dimensions to see that on the horizon 

g ab VJ V J = 0. (4.9) 

That is, the 4-vector V/ is null on the horizon. In fact, on the horizon, where in terms 
of the 3-normal n we can decompose v# = c s n + vm, we can furthermore write 

l|V/|| 



(V/r„ = (g ab V b f) 



ii 



tt 2 H c H 



(i;v„) 



H 



(4.10) 



That is, not only is the 4-vector V/ null on the horizon, it is also a 4-tangent to 
the horizon — so (as in general relativity) the horizon is ruled by a set of null curves. 
Furthermore, extending the 3-normal n to a region surrounding the horizon (for instance 
by taking n = V//|| V/||) we can quite generally write v = v± n + vn. Then away from 
the horizon 



9 



n 2 d 



(4.11) 



That is, the 4-vector V/ is spacelike outside the horizon, null on the horizon, and 
timelike inside the horizon. 

Stationary but non-static Killing horizons: If the stationary horizon is Killing, 
then even if we do not explicitly know what the relevant Killing vector K a is, we know 
that its norm has to vanish on the horizon, and so the norm of this horizon-generating 
Killing vector is of the form 

„,2 



K a K a = Q (cj 



-Qc 2 s 



(4.12) 



for some unknown (but for current purposes irrelevant) function Q. Following closely 
the argument for the static case, mutatis mutandis, we have [3-5] 



9h 



n ■ V(c 2 - v\ 
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n ■ Vfc, 



H 



d{c s 



dn 



(4.13) 



H 



which can also be compactly written in terms of the horizon-crossing Mach number, 
Mj_ = v±/c s , as 













dM ± 


9h = c 2 H 


n • V(v ± /c s ) 




n • VM ± 




dn 



(4.14) 
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We again emphasise that this already works in the geometric acoustics framework, and 
that there is no need to make the more restrictive assumptions corresponding to wave 
acoustics that were made in references [3] and [5]. If the horizon is non-Killing then 
one must be even more devious. 

Stationary but non-static non-Killing horizons: Such situations are, from a 
technical perspective, much more difficult to deal with. Such behaviour cannot occur 
in standard general relativity, where the Einstein equations stringently constrain the 
allowable spacetimes, but there seems no good reason to exclude it for acoustic horizons. 
Unfortunately, when it comes to explicit computations of the surface gravity there are 
still some unresolved technical issues. There is still a lot of opportunity for significant 
new physics hiding in these non-Killing horizons. 

5 Relativist ic acoustics 

Full relativistic acoustics (either special relativistic or general relativistic) adds a few 
other quirks which I briefly describe below. (See early astrophysical work by Mon- 
crief [26], a more recent cosmological framework developed in [27], and a pedagogical 
exposition in reference [7] for details.) Note that the interest in, and need for, rel- 
ativistic acoustics is driven by astrophysical and cosmological considerations, not by 
direct laboratory applications. There are at least three situations in which relativistic 
acoustics is important: 

• Speed of sound comparable to that of light. 

In any ideal gas once kT 3> m c 2 then p w |p and so c s ~ -^c. 

This is physically relevant, for instance, in various stages of big bang cosmology. 

• Speed of fluid flow comparable to that of light. 

This is physically relevant, for instance, in some black hole accretion disks and/or 
the jets emerging from active galactic nuclei (AGNs). 

• Tight binding: p < g or <C moc 2 . 

Once the pressure is an appreciable fraction of the energy density, or the absolute 
value of the chemical potential is much smaller than the rest mass, then the usual 
derivation of the conformal factor appearing in the wave version of the acoustic 
metric must be significantly modified. 

This is physically relevant, for instance, in cores of neutron stars. 

It is somewhat unclear at present as to whether relativistic acoustics can be made 
directly relevant for laboratory physics. Some first steps in this regard may be found 
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in reference [28] , where the possibility of experimentally constructing relativistic BECs 
is considered. 



Geometric acoustics: If one works with special relativistic acoustics, rather than 
non-relativistic acoustics, then at the level of ray acoustics one will simply obtain an 
acoustic variant of the Gordon optical metric 



The only difference is that the refractive index has now been replaced by the ratio of 
the speed of sound to the speed of light: n~ l (x) —> c s (x)/c. (The 4- velocity of the 
medium is still V a (x), and the conformal factor Q(x) is still undetermined.) In general 
relativistic acoustics this would become 



where f a b{x) is now the general relativistic physical background metric obtained by 
solving the Einstein equations, and g ab is the acoustic metric for the acoustic pertur- 
bations in the fluid flow. Note that the 4- velocity V a (x) now has to be a timelike unit 
vector with respect to the background metric f a b{x). For ray acoustics this is all one 
can say. 

Wave acoustics: One can again go to wave acoustics, deriving a wave equation by 
linearizing the general-relativistic version of the Euler equations. The same sort of 
technical assumptions must be made, (irrotational, barotropic, and inviscid), and one 
now obtains a slightly more complicated formula for the conformal factor [7] 



Here n is the number density of particles, and g is the energy density (rather than the 
mass density p), while c s is the speed of sound. The quantity p is the pressure, and d 
is the number of space dimensions. 

Non-relativistic limit: In the non-relativistic limit among other things we certainly 
have p <C g. Also in terms of the average particle mass fn one has 



g ab = n 2 (rj ab + V a V b ) - -§ VaVb 







(5.3) 



Q 




(5.4) 



and so 



n 




(5.5) 



c s (g + p) 



c s g 



c s (mc 
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thereby (as required for internal consistency) reproducing the correct limit for the 
conformal factor. 

The correct limit for the tensor structure is more subtle. (A suitable discussion 
can be found in reference [7].) Formally taking the limit c — > oo, but holding c s and v 
fixed, a brief calculation yields: 



9oo = ^ 2 



■1 + 7 



goi = -tt 2 
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c 2 
cr 



-7 2 

2 ' 



-Si 2 ^ 

CT 



+ ... 
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g oi = fi 2 <5. 



(5.6) 

(5.7) 
(5.8) 



Then, switching from (ct, x) coordinates to (t, x) coordinates, the relativistic g a b of 
this section correctly limits to the non-relativistic g a b of the previous section. 



6 Bose— Einstein condensates 

Bose— Einstein condensates (BECs) provide a particularly interesting analogue model 
because they are relatively easy to construct and manipulate in the laboratory, and 
specifically because the speed of sound is as low as a few centimetres per second. Most 
work along these lines has focussed on non-relativistic BECs. Suitable background 
references are [29-35]. See also the companion chapter by Balbinot et al. in the current 
volume [36]. In view of the coverage of this topic already provided in that chapter, I 
shall not have more to say about it here. 

In contrast, I will briefly discuss the relativistic BEC model of Fagnocchi et al. 
that is presented in reference [28]. While relativistic BECs do not seem currently to 
be a realistic experimental possibility, the theoretical treatment introduces some new 
issues and effects. The relativistic BECs naturally lead to two quasiparticle excitations, 
one massless and one massive, with rather complicated excitation spectra. (In this 
sense the relativistic BECs are reminiscent of the "massive phonon" models obtained 
from multiple mutually interacting non-relativistic BECs [37-39].) In the relativistic 
BEC one obtains a 4th-order differential wave equation for the excitations, which is 
ultimately why one has two branches of quasiparticle excitations. In the limit where 
the relativistic generalization of the so-called quantum potential can be neglected, the 
wave equation simplifies to the d'Alembertian equation — for a relativistic acoustic 
metric of the generalized Gordon form discussed in the previous section. In the limit 
where both relativistic effects and the quantum potential can be neglected, one recovers 
the (wave acoustic version of) Unruh's non-relativistic acoustic metric. 
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Madelung representation: There is an important and non-obvious technical point 
to be made regarding the linearization of the Madelung representation in a BEC con- 
text, (or in fact in any situation where one is daeling with a non-linear Schrodinger-like 
equation). For any complex field ip the Madelung representation is 

V = Vpe*. (6.1) 

When linearizing, (which is the basis of separating the system into background plus 
excitation, or condensate plus quasiparticle), there are at least three things one might 
envisage doing: 

1. Take ip = ip + e fa + 0(e 2 ). 

2. Take p = p + e p x + 0(e 2 ), and = O + e 0i + 0(e 2 ). 

3. Take tfj = ^ {1 + e X + 0(e 2 )}. 
Note that routes 2 and 3 are related by: 

Pi X + X ] , x-x ] 



po 2 ' r 2i 



(6.2) 



Mathematically, all three routes must carry the same intrinsic physical information, 
but the clarity with which the information can be extracted varies widely depending 
on the manner in which the perturbative analysis is presented. When actually carrying 
out the linearization, it turns out that route 1 is never particularly useful, and that 
routes 2 and 3 are essentially equivalent for a non-relativistic BEC, ultimately leading 
to formally identical wave equations. In contrast, for relativistic BECs it is route 3 
that leads to the cleanest representation [28], while route 2 leads to a bit of a mess [3]. 
(A mess involving integro-differential equations.) This is not supposed to be obvious, 
and will not be obvious unless one tries to carefully work through the relevant technical 
literature. With hindsight, route 3 appears to be the superior way of organizing the 
perturbative calculation. 



7 Surface waves and blocking horizons 

Surface waves (water-air, or more generally waves on any fluid-fluid interface) are de- 
scribed by an incredibly complex and subtle theoretical framework — one of the major 
technical complications comes from the fact that surface waves are highly dispersive, 
with a propagation speed that is very strongly frequency dependent. Thus, insofar as 
one can put surface wave propagation into a Lorentzian metric framework, one will 
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have to adopt a "rainbow metric" formalism with a frequency dependent metric. The 
trade-off is that this system is relatively easily amenable to laboratory investigation 
through "wave tank" technology [40-42]. 

Surface waves in the geometric limit: As long as the wavelength and period of 
the surface wave are small compared to the distances and timescales on which the depth 
of water is changing one can usefully work in the geometric (ray) limit. Under those 
conditions one can write a 2+1 dimensional metric to describe ray propagation: 



ds 2 



n 2 {- C 2 SW dt 2 + (dx - v dty} = n 2 



{-(cl->- 



2v-dxdt+ ||dx|| 2 }. (7.1) 



Here c sw is the speed of the surface waves in the comoving frame (that is, comoving with 
the surface of the fluid), and v is the (horizontal) velocity of the surface. Unfortunately 
the speed c sw is a relatively complicated function of (comoving) frequency, depth of the 
water, density of the fluid, the acceleration due to gravity, the surface tension, etcetera. 
(See for instance references [43-46].) 

Now in terms of a 3 x 3 matrix, this is equivalent to defining the metric tensor 
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The indices a, b, c, take on values in {0,1,2}, corresponding to both time and 
(horizontal) space, whereas indices such as i, j, k, . . . , take on values in {1, 2}, corre- 
sponding to (horizontal) space only. The corresponding inverse metric is 
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(7.3) 



In the fluid dynamics community, one most often restricts attention to 1 spatial di- 
mension, then surface waves are said to be "blocked" whenever one has ||v|| > c sw , 
and one will encounter considerable attention paid to this concept of "wave blocking" 
in that community. This is what a general relativist would instead call "trapping", 
and consequently the mixed terminology phrase "blocking horizon" has now come into 
use within the analogue spacetime community. Note that instead of speaking of Mach 
number (appropriate for acoustic propagation through the bulk of a medium), in a sur- 
face wave context it is the Froude number that governs the formation of ergo-regions 
and horizons. 

Surface waves in the physical limit: Moving beyond the geometric/ray approx- 
imation for surface waves is mathematically rather tricky. Within the fluid dynamics 
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community relevant work is based on the Boussinesq approximation [47, 48], and its 
modern variants [49]. Within the analogue spacetime community, see particularly the 
basic theoretical work in reference [50], and in the related chapter [51] in this volume. 
(See also [52, 53] for a more applied perspective.) Physically, in addition to the presence 
of dispersion, a second complicating issue is this: The fluid at the surface is moving 
both vertically (the wave) and horizontally (the background flow), while at the base 
of the fluid (which may be at variable depth), the no-slip boundary condition enforces 
zero velocity. 

Based on the three-dimensional Euler and continuity equations one then has to 
construct an interpolating model for the fluid flow that connects the surface motion to 
the zero-velocity motion at the (variable depth) base. Once this is achieved, one throws 
away the interpolating model and concentrates only on the physical observable: the 
motion of the surface. The analysis is mathematically and physically subtle, and (in the 
physical or wave limit) the theoretical framework for surface waves is nowhere near as 
clean and straightforward as for barotropic inviscid irrotational acoustic perturbations 
travelling through the bulk. 

Experiments: The key benefit of surface waves is that the propagation speed c sw 
is easily controllable by adjusting the depth of fluid, that background flows are easily 
set up by simple mechanical pumps, and that "wave tank" and related technologies 
are well understood and well developed. (See for example, the early 1983 experiments 
by Badulin et al. [40].) This particular analogue spacetime has recently led to several 
very interesting experimental efforts [41, 42, 54]. For instance, Weinfurtner et al. have 
performed an experiment looking at the classical (stimulated) analogue of Hawking 
radiation from a blocking horizon, and have detected an approximately Boltzmann 
spectrum of Hawking-like modes [41], while Rousseaux et al. have experimentally in- 
vestigated the related "negative-norm modes" [42] . The relation between the "hydraulic 
jump" and blocking horizons has been experimentally investigated by Jannes et al. [54]. 
Some related theoretical developments are reported in [45, 55]. Work on this topic is 
ongoing. 

8 Optical fibres/optical glass 

In an optical context, related "optical blocking" phenomena occur when a "refractive 
index pulse" (RIP) is initiated in an optical fibre [56], or in optical glass [57-59]. The 
basic idea is that things are arranged so that while the RIP moves at some speed 
vrip, the velocity of light outside the RIP is greater than the velocity of the RIP 
c outside > ^rip, while inside the RIP we have the contrary situation c insi de < v-^ip. This, 
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(certainly within the geometric optics framework), sets up a "black" horizon at the 
leading edge of the RIP, and a "white" horizon at the trailing edge. (For technical 
details see references [56-59].) Some subtleties of the theoretical analysis lie in the 
distinction between group and phase velocities — are we dealing with "phase velocity 
horizons" or "group velocity horizons" ? Other technical subtleties have to do with the 
transition from geometric optics to wave optics — there are a number of complex and 
messy technical details involved in this step. 

An intriguing experiment based on these ideas has been carried out by Belgiorno et 
ai, with results reported in reference [60]. While it is clear that some form of quan- 
tum radiation has been detected, there is still some disagreement as to whether this is 
(analogue) Hawking radiation, or possibly some other form of quantum vacuum radia- 
tion [61-63]. Work on this topic is ongoing. 

9 Other models 

A complete and exhaustive catalogue of other analogue models would be impractical. 
See the Living Review article on "Analogue gravity" for more details [3]. Selected 
models, (a necessarily incomplete list), include: 

• Electromagnetic wave guides [64]. 

• Graphene [65, 66]. 

• Slow light [67-71]. 

• Liquid helium [72, 73]. 

• Fermi gasses [74, 75]. 

• Ion rings [76]. 

Beyond the issue of simply developing analogue models, there is the whole subject of 
using analogue models to probe, (either theoretically or more boldly experimentally), 
a whole raft of physics questions such as directly verifying the existence of Hawking 
radiation, the possibility of Lorentz symmetry violations [77], the nature of the quantum 
vacuum [78-80], etcetera. For more details, see reference [3], and other chapters in this 
volume. 
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10 Discussion 



The general theme to be extracted from these considerations is this: The propagation 
of excitations (either particles or waves) over a background can often (not always) be 
given a geometric interpretation in therms of some "analogue spacetime" . As such a ge- 
ometric interpretation exists, there is a strong likelihood of significant cross-fertilization 
of ideas and techniques between general relativity and other branches of physics. Such 
possibilities have increasingly attracted attention over the last decade, for many rea- 
sons. The other chapters in these proceedings will explore these ideas in more specific 
detail. 
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